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Abstract 

Plant cells are routinely exposed to various pathogens and environmental stresses that cause cell wall perturbations. 
Little is known of the mechanisms that plant cells use to sense these disturbances and transduce corresponding 
signals to regulate cellular responses to maintain cell wall integrity. Previous studies in rice have shown that 
removal of the cell wall leads to substantial chromatin reorganization and histone modification changes 
concomitant with cell wall re-synthesis. But the genes and proteins that regulate these cellular responses are still 
largely unknown. Here we present an examination of the nuclear proteome differential expression in response to 
removal of the cell wall in rice suspension cells using multiple nuclear proteome extraction methods. A total of 382 
nuclear proteins were identified with two or more peptides, including 26 transcription factors. Upon removal of the 
cell wall, 142 nuclear proteins were up regulated and 1 12 were down regulated. The differentially expressed 
proteins included transcription factors, histones, histone domain containing proteins, and histone modification 
enzymes. Gene ontology analysis of the differentially expressed proteins indicates that chromatin & nucleosome 
assembly, protein-DNA complex assembly, and DNA packaging are tightly associated with cell wall removal. Our 
results indicate that removal of the cell wall imposes a tremendous challenge to the cells. Consequently, plant cells 
respond to the removal of the cell wall in the nucleus at every level of the regulatory hierarchy. 
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Background 

The cell wall is a critical extracellular structure that pro- 
vides protection and structural support in plant cells. It 
controls the cell shape and allows the turgor pressure to 
build up and maintain an upright position for plants. In 
addition, it glues the cell together and serves as a barrier 
for pathogen infection and insect and animal damage. 
Plant cells are routinely exposed to various pathogens 
and environmental stresses that cause cell wall perturba- 
tions. Insect and herbivore bites and wind are common 
factors contributing to cell wall damage. Little is known 
about the mechanisms that plants use to sense these 
disturbances and transduce the signals to stimulate 
responses to maintain cell wall integrity. It has been 
demonstrated in yeast cells that transient damage to cell 
wall leads to induction of cell wall-related genes as a 
compensatory response to maintain cell integrity [1]. 
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However, in spite of clues from many stress -related stud- 
ies, it is unknown if such a mechanism exists in plant 
cells. 

Plant cells can rapidly re-synthesize the cell wall after 
the cell wall is removed [2], The plant protoplast culture 
is an excellent experiment displaying the astonishing cell 
wall re-synthesis capability. Interestingly, the cell wall 
re-synthesis mechanism in protoplasts is probably differ- 
ent from the one used for new cell wall synthesis during 
cell division [3]. Tan et al. (2011) found that removal of 
cell wall leads to cell wall synthesis at multiple sites in 
protoplasts [3]. In contrast, new cell wall synthesis dur- 
ing cell division is limited to only one site- the cell plate 
derived from the phragmoplast. In addition, substantial 
chromatin reorganization was observed in protoplasts. 
The chromatin reorganization was associated with his- 
tone modification changes at multiple modification sites 
of histones as shown in Western blot studies with mul- 
tiple histone modification specific antibodies. The his- 
tone acetylation changes at H3K18 and H3K23 following 
cell wall removal and regeneration were further verified 
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and quantified using isotope labeling assisted mass spec- 
trometry analysis. In addition, 136 up regulated and 94 
down regulated proteins were identified using shot gun 
proteomics and label-free quantification analysis [3]. 
Sharma et al (2011) examined the transcriptome re- 
sponse to enzymatic removal of cell wall [4]. They found 
that kinases, transcription factors and genes predicted to 
be involved in cell wall-related functions were enriched 
in the differentially regulated gene category. In addition, 
rice lines carrying Tosl7 mutations in genes up- 
regulated during cell wall removal exhibited dwarf phe- 
notypes. Many of the genes up-regulated during cell wall 
re-synthesis following cell wall removal are also up- 
regulated in response to infection and environmental 
perturbations, indicating that there is a coordinated 
response to diverse types of stress. 

The nucleus is the most prominent organelle that con- 
tains majority of the genetic materials in eukaryotes. It is 
the site of DNA replication, RNA transcription, and 
ribosome preassembling. The nucleus is surrounded by a 
double membrane called the nuclear envelope. The nu- 
cleus contains several subcompartments [5], including 
nucleolus, euchromatin domains, heterochromatin do- 
mains, cajal bodies, speckles, and other domains. The 
nuclear matrix is a karyoskeletal, non-histone structure 
that serves as a support for the genome and nuclear ac- 
tivities [6]. The essential roles of the nuclear activities to 
the cell suggest that the nucleus is the most important 
control center of the cell. 

The nuclear proteome is highly complicated, with pro- 
teins ranging from very low copy transcription factors to 
highly abundant core histone proteins and ribosomal 
proteins. In plants, the nuclear proteome has been ex- 
amined by several laboratories in different organisms. 
The nuclear proteins were extracted using different 
methods for proteomics studies, including Trizol extrac- 
tion [7], fractionation with differential ionic strength [8], 
high NaCl concentration [9], HEPES buffer [10], lysis 
buffer [11,12], and phenol extraction [13]. In rice, 
glucose-responsive nuclear proteins were extensively 
examined [9]. Nuclear enriched proteomes were also 
studied in different tissues in rice [9,12,13]. The nuclear 
proteome response to cold stress has been examined in 
Arabidopsis with several transcription factors shown to 
be differentially regulated under stress. Nucleolar, nu- 
clear matrix, and nuclear pore complex proteomes were 
also examined in Arabidopsis [14-17]. Although many 
nuclear proteome studies have been reported, the num- 
ber of low abundance transcription factors identified in 
each study was usually less than ten. When nuclei- 
enrichment was combined with a DNA binding affinity 
column, about a dozen transcription factors were identi- 
fied [9], suggesting that improving the nuclear protein 
purification and extraction methods may lead to a better 



coverage of the nuclear proteome, particularly the low 
abundance proteins. 

Although differential histone modifications and chro- 
matin reorganization in response to cell wall removal 
and regeneration have been observed in rice, the regula- 
tory network controlling the process is still largely un- 
known. No regulatory genes specifically involved in this 
process have been identified at the protein level. In this 
report, we used multiple nuclear proteome extraction 
methods to examine the nuclear proteome response to 
the removal of the cell wall. A large number of nuclear 
proteins including histone modification proteins, chro- 
matin structure regulatory proteins, and transcription 
factor proteins were identified. Our studies substantially 
advanced our understanding of the plant nuclear prote- 
ome and cellular responses to cell wall removal. 

Results 

Cell wall removal stimulates active cell wall synthesis 

To study how plant cells respond to the disturbance of 
cell wall, we examined cellular responses to the enzym- 
atic removal of cell wall using rice suspension culture 
cells, the OC cell line [18,19]. Because of the unique cell 
wall structure of plants in the grass family, multiple 
hours of enzyme digestion are required to completely re- 
move the rice cell wall [20-22]. After 9 hours of enzyme 
digestion, the cell wall was completely removed as re- 
vealed by the stain with Fluorescent Brightener 28, a fluor- 
escent dye with specific polysaccharide binding activities 
(Figure 1A). Followed by 2 to 4 hours of culture of the 
protoplasts, new cell wall started to emerge (Figure IB). 
After 48 hours of culture, the relatively spherical and 
smooth surface of protoplasts changed, suggesting the 
recovery of cell wall, which is verified by Fluorescent 
Brightener 28 stain (Figure 1C). We found that over 90% 
protoplasts could regrow their cell wall, suggesting that 
our protoplast isolation and culture is an excellent system 
to examine cellular response to the removal of cell wall. 

Nuclei enrichment and assessment 

A high quality and large scale purification of nuclei is 
vital to nuclear subproteome analysis. We obtained nu- 
clei from protoplasts and suspension cells, respectively. 
DAPI-staining of purified nuclei fractions from both sus- 
pension cells and protoplasts revealed that we obtained 
nuclei in a large scale from both suspension cells and 
protoplasts without any clear contamination from organ- 
elles such as chloroplasts and mitochondria as observed 
under the microscope (Figure 2A & B and data not 
shown). We validated the nuclear enrichment by Western 
blots with antibodies specific for known nuclear and cyto- 
solic proteins. The Western blot results showed that 
histone H4 was highly enriched in the nuclear fraction 
compared to total protein extraction when equal amount 
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Figure 1 Microscopy images of cultured rice protoplasts (from suspension cells) following the cell wall regeneration time course. CLSM 
was used to observe the protoplasts stained by a fluorescent dye, Fluorescent Brightener 28, with polysaccharide specific binding activities. The 
excitation wavelength at 492 nm and emission at 520 nm were used. The protoplast culture times are Ohrs (A), 4 hrs (B), and 48 hrs (C). The 
arrows point at the positions of cell wall syntheses. The magnification is revealed by the scale bar. 



of proteins were loaded (Figure 2C). In contrast, cytosolic 
fructose- 1, 6-bisphosphatase (cFBPase) and vacuolar 
protein VHA-E were only detected in the total protein 
fraction (Figure 2C), indicating that the nuclear proteins 
were successfully enriched. 

Comparison of nuclear protein extraction methods 

Nuclear subproteomes have been studied with different 
protein extraction methods, including Trizol extraction 
[7], fractionation with differential ionic strength [8], high 
NaCl concentration [9], HEPES buffer [10], lysis buffer 
[11,12], and phenol extraction [13]. However, the low 
abundant nuclear proteins identified by mass spectrom- 
etry are still limited in plants. To optimize the method 
for nuclear protein identification, we tested different 
nuclear proteome extraction and fractionation methods 
as revealed in Figure 3A. To determine if a protein was 
localized in the nucleus, GO annotations were obtained 
from GORetriever, a tool available at AgBase [23,24]. 
We found that a combination of the phenol extraction 
with acid re-extraction could improve the nuclear sub- 
proteome coverage (Figure 3B). Phenol extraction of the 
nuclei derived from protoplasts and suspension cells 
followed by LC-MS/MS identified 251 and 115 nuclear 
proteins, respectively. Acid extraction followed by LC-MS 
/MS identified 137 and 165 nuclear proteins, respectively. 
When the phenol extracted samples were re-extracted by 
sulfuric acid and examined with LC-MS/MS, 113 and 144 
nuclear proteins were identified in the nuclear samples of 
protoplasts and suspension cells, respectively. Among 
them, 15 and 47 proteins, respectively, were new proteins 
that were not identified by either the phenol or acid 
extraction method. The total nuclear proteins identified by 
each of the extraction methods are listed in Additional 
file 1: Table SI. Due to some overlap, overall we identified 
382 nuclear proteins with two or more peptides. Among 



them, 26 were transcription factors. All proteins discussed 
and presented in this study met the criterion of two or 
more matched peptides. To verify our protein identifica- 
tion results, a reverse database of O. sativa was searched 
using the reverse database functionality in Bioworks 3.2 as 
previously reported (Additional file 1: Table S2) [3]. The 
peptide false discovery rate (FDR) for the entire dataset 
was 0.58%, while the protein FDR was 1.51%. 

Analysis of the total identified peptides showed that 
about 31% of the peptides identified using phenol extrac- 
tion were nuclear protein peptides. When the sample 
was re-extracted by acid, 67% of the identified peptides 
were nuclear protein peptides. Nine of the top 10 most 
abundant proteins (based on peptide counts) identified 
in the acid re-extraction samples were histones (Table 1). 
In contrast, none of the 10 most abundant proteins 
extracted by phenol alone were histones although the 
majority was nuclear proteins (Table 1), suggesting that 
acid re-extraction enriched nucleic acid associated pro- 
teins. Meanwhile, 47% of the peptides identified in sam- 
ples directly extracted by acid were nuclear protein 
peptides. Of the 10 most abundant proteins identified by 
acid extraction, three were histones and three were 
nucleolar proteins. 

Differentially expressed proteins in response to cell 
wall removal 

Upon removal of cell wall, rice cells display substantial 
chromatin decondensation and reorganization [3]. To 
identify nuclear proteins that may be involved in chro- 
matin decondensation and reorganization; we examined 
differentially expressed nuclear proteins upon the re- 
moval of cell wall. To reveal the differentially expressed 
proteins, we compared the suspension cell nuclear prote- 
ome with the protoplast nuclear proteome extracted by 
phenol extraction, acid re-extraction, and acid extraction, 
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Figure 2 Microscopy images of isolated rice (O. sativa) suspension cell and protoplast nuclei and Western blot analysis of purified nuclear 
proteins. (A) Image of purified suspension cell nuclei after DAPI staining. (B) Image of protoplast with cell wall regeneration (4 hrs) nuclei after DAPI 
staining. A small volume of the purified nuclei was stained with DAPI (0.5 ug/ml) for 5 minutes and images were taken under a DAPI-filter. The 
magnification is revealed by the scale bar. (C) Nuclei enrichment revealed by Western blots. Antibodies against H4, V-ATPase, E, and cFBPase were used to 
assess the protein quantity in the total protein fraction and suspension cell and protoplast nuclei, respectively. 20 ug of proteins were loaded in each lane. 
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Figure 3 Protein extraction methods utilized in this study and nuclear proteins identified in the various fractions. (A) Suspension cell 
nuclei and protoplast nuclei were extracted with phenol alone, phenol and 0.4 N sulfuric acid (A & P), and 0.4 N sulfuric acid alone, respectively, 
followed by subsequent mass analysis. (B) Color-scheme Venn diagrams revealing identified nuclear proteins in each extraction procedure and 
the overlap among extraction procedures in suspension cell and protoplast nuclear samples. The numbers in circle areas equal the protein 
number identified. Purple: Acid extraction; Green: Phenol Extraction; Blue: Phenol-Acid double extraction. 



respectively. A non-labeling quantification method was 
used for differential regulation analysis. Previous reports 
and our studies have shown that the spectral count and 
ZX corr score methods generated identical results in all 
studies [3,25]. But the SX corr score method provided 
values for direct comparison of protein fold-change. There- 
fore we used the SXc 0rr score method. Also, the sum of 
SEQUEST ZXcorr has been shown to compare suitably with 
the concentrations of a known protein mixture in serial di- 
lutions [26]. In the ZXcorr method a preliminary list is built 
using all scans for peptides with an Xcorr (generated by 
TurboSEQUEST™ (Bioworks Browser 3.2, Thermo Elec- 
tron)) above the threshold used for protein identification. 
Finally, the values of suspension cell control ZXcorr versus 
protoplast treatment ZX cor r for individual proteins identi- 
fied are compared and statistically significant changes are 
used to assign regulation and fold-change [25,26]. The 
Xcorr values generated from TurboSEQUEST were used for 
ZXcorr quantification as reported by Nanduri and Bridges 
[26,27], in which three biological replicas of each sample 
treatment is required. The quantitative analysis criteria and 
procedure were identical to previously reported [3,25,27]. 
Differential expression was only considered for proteins 
with a p-value < 0.05. 



Following removal of cell wall, 142 nuclear proteins 
with a p-value < 0.05 displayed differential up regulation 
and 112 nuclear proteins with a p-value < 0.05 displayed 
differential down regulation (Additional file 1: Table S3). 
To validate the protein differential expression results gen- 
erated by the ZX cor r method between the suspension cells 
and protoplast at the transcriptional level, we randomly 
selected nine differentially expressed proteins for RT-PCR 
and real-time PCR analysis (Figures 4 and 5). The expres- 
sion levels of these genes correlated with the non-labeled 
protein quantification results, providing further support 
for our protein quantification results. To further analyze 
the differentially regulated proteins, functional classifica- 
tion of the differentially expressed nuclear proteins was 
carried out according to the gene ontology (GO) rules 
using AgBase at http://www.agbase.msstate.edu/ [24] and 
ortholog and Pfam domain information available for all 
proteins identified with two or more peptides was col- 
lected using the tools provided by the TIGR Rice Genome 
Annotation Project (http://rice.plantbiology.msu.edu/). 
Ortholog and Pfam domain information available for the 
identified proteins is presented in Additional files 2 and 3, 
respectively. Three independent gene ontologies were 
used to describe the function of gene products such as 
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Table 1 The most abundant proteins identified in phenol, acid, and phenol-acid extracted suspension cell and 



protoplast nuclear samples 

TIGR ID Annotation Peptide (Hits) 

Phenol-Extracted 

LOC_Os03g22740 Nucleolar protein N0P5-1, putative, expressed 108 

LOC_Os08g 04240 Cysteine-rich repeat secretory protein 55 precursor, putative, expressed 75 

L0C_0s1 1 g 1 0480 Dehydrogenase, putative, expressed 73 

LOC_Os03g22730 Nucleolar protein NOP5-1, putative, expressed 67 

LOC_Os05g08360 rRNA 2-O-methyltransferase fibrillarin 2, putative, expressed 51 

LOC_Os04g40950 Glyceraldehyde-3-phosphate dehydrogenase, putative, expressed 48 

LOC_Os08g 04250 Cysteine-rich repeat secretory protein 55 precursor, putative, expressed 47 

LOC_Os03g22880 Nucleolar protein 5A, putative, expressed 46 

LOC_Os02g57590 rRNA 2-O-methyltransferase fibrillarin 2, putative, expressed 44 

LOC_Os02g38920 Glyceraldehyde-3-phosphate dehydrogenase, putative, expressed 44 
Acid-Extracted 

LOC_Os08g 04240 Cysteine-rich repeat secretory protein 55 precursor, putative, expressed 125 

LOC_Os08g 04250 Cysteine-rich repeat secretory protein 55 precursor, putative, expressed 107 

LOC_Os08g04210 Cysteine-rich repeat secretory protein 55 precursor, putative, expressed 81 

LOC_Os04g52960 Nucleolin, putative, expressed 70 

LOC_Os03g22730 Nucleolar protein NOP5-1, putative, expressed 65 

LOC_Os03g22740 Nucleolar protein NOP5-1, putative, expressed 63 

LOC_Os07g44190 h/ACA ribonucleoprotein complex subunit 4, putative, expressed 51 

LOC_Os01g61920 Histone H4 42 

LOC_Os05g38640 Probable histone H2A.4 41 

LOC_Os05g02300 Probable histone H2A.6 41 
Phenol-Acid-Double Extracted 

LOC_Os02g56960 Ribosomal protein, putative, expressed 146 

LOC_Os05g38640 Probable histone H2A.4 99 

LOC_Os05g02300 Probable histone H2A.6 99 

LOC_Os03g 17100 Probable histone H2A.5 99 

LOC_Os07g36500 Histone H4 70 

LOC_Os01g61920 Histone H4 70 

LOC_Os01g05900 Histone H2B.10 43 

LOC_Os01g05630 Histone H2B.4 43 

LOC_Os01g05610 Histone H2B.3 43 

LOC_Os07g36140 Probable histone H2A.2 39 



cellular component (CC), molecular function (MF) and 
biological process (BP) [28]. GO annotations were obtained 
from GORetriever, a tool available at AgBase [23,24], GO 
classification was carried out using tools available at 
AgBase [23,24] and AgriGO [29,30]. Functional classifica- 
tion for differentially regulated proteins in categories as 
cellular component, molecular function, and biological 
process were found. The results are presented in Figure 6 
and Additional file 1: Figure SI A and SIB, respectively. 

Figure 6 shows the significantly enriched GO biological 
processes of differentially expressed nuclear proteins. The 



biological processes tightly associated with cell wall regener- 
ation included chromatin assembly, nucleosome assembly, 
macromolecular complex subunit organization, protein- 
DNA complex assembly, and DNA packaging (Figure 6). 

Differential expression of transcriptional regulation proteins 

Identifying regulatory proteins such as transcription 
factors controlling cellular response to cell wall removal 
is essential for revealing the cellular regulatory network. 
However, transcription factors are difficult to detect 
by mass spectrometry due to low copy numbers. We 
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Figure 4 RT-PCR analysis of differentially expressed genes in suspension cell and protoplast nuclei. Rice ubiquitin gene was used as an 
internal control. Equal amount of cDNA template was used for suspension cell and protoplast nuclei cDNA samples. 35 cycles were used. Primers 
utilized in this study are provided in Additional file 4. 



successfully identified 26 transcription factors and found 
that several of them were differentially regulated, including 
multiple zinc finger proteins. While zinc finger C-x8-C- 
x5-C-x3-H type family proteins LOC_Os02g06584 and 
LOC_Os06g46890 proteins were up regulated in response 
to cell wall removal, zinc finger family protein (LOC_Os04 
g57010) and ZOS3-23-C2H2 zinc finger protein (LOC_0 
s03g61640) were down regulated. Other differentially reg- 
ulated transcription factors included Whirly transcription 



factor domain containing protein (LOC_Os06g05350), 
Helix-loop-helix DNA binding domain containing protein 
(LOC_Os02g39140), transcription factor TF2 (LOC_Os05 
g03740), and putative transcription factor (LOC_Os09g 
27850). Other proteins that might be involved in transcrip- 
tional regulation were also differentially regulated. The 
SKIP (SKI-interacting protein) is an essential spliceosomal 
component and transcriptional co-regulator, which may 
provide regulation by coupling transcription initiation with 
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Figure 5 Quantitative real-time PCR analyses of differentially expressed genes in suspension cell and protoplast nuclei. Equivalent 
amount of cDNA template was used for each sample and rice ubiquitin qene was used as an internal control. 
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Figure 6 Enriched GO biological processes of differentially expressed nuclear proteins. Figure displays the significantly enriched biological 
processes revealed by GO annotation analysis for differentially expressed proteins. The top line in each box is the GO identifier of the term and 
statistical significance (multiple hypothesis corrected p-value, lower is more significant) of that annotation. The middle line in each box is a 
description of the GO term. The four numbers on the bottom line are the number of nuclear proteins that had this annotation, the number of 
nuclear proteins that had any annotation (192), the total number of proteins that had the annotation, and the total number of proteins that had 
any annotation (24460). The color of the box indicates the significance of the term as indicated by the legend on the bottom left corner. White 
boxes are not significant. 



splicing [31]. A SKIP/SNW domain containing protein 
(LOC_Os02g52250) was down regulated. The BRCA1 C 
terminus domain containing protein (LOC_Os03g49210) 
was also up regulated. SSRPl-like FACT complex subunit 
(LOC_Os01g08970) was found to be up regulated. The 
FACT complex contains proteins such as SSRP1 and 
Sptl6, which are connected with transcriptional elongation 
[32,33]. Finally, a putative DNA-directed RNA polymerase 
subunit (LOC_Os09g02284) was also up regulated. RNA 
polymerase II is a multi-subunit holoenzyme composed of 
ten to twelve protein subunits (RPB1-RPB12) [34-36]. 
LOC_Os09g02284 is an ortholog of the Arabidopsis pro- 
tein DNA-directed RNA polymerase II subunit RPB3-B 
(At2gl5400), which composes the core element of the 
RNA polymerase II protein. 

Differential expression of chromatin structure and 
modification proteins 

A large number of genes regulating chromatin structure 
and function were differentially regulated, including core 
histone proteins, core histone domain containing pro- 
teins, HMG proteins, histone modification proteins, and 
nucleosome remodeling proteins. Interestingly, several 
core histone domain containing proteins were up regu- 
lated. The function of these genes remains to be further 
explored. Meanwhile, the H3 proteins were also up regu- 
lated. The HMG proteins are present in all tissues of 
eukaryotes, leading many to believe that HMG proteins 
are central for proper cellular function [37]. Johns et al. 
(1982) [38] estimated that HMG proteins bind to <10% 
of the nucleosomes, making them the second most 
abundant family of chromosomal proteins with probable 



structural function in the nucleus [37]. While HMG-Y- 
related protein (LOC_Os09g23730) was up regulated, 
the putative HMG1/2 (LOC_Os06g51220) was down 
regulated. 

The histone modification proteins are believed to 
regulate the access of transcription factors, chromatin 
modifying enzymes, and chromatin remodeling factors 
to nucleosomal DNA by chemical modifications to his- 
tones. We found that while the putative histone deacetylase 
(LOC_Os07g06980) was up regulated, the histone methyl- 
transferase (LOC_Os05g41172) was down regulated. A 
histone lysine N-methyltransferase H3 lysine-9 specific 
SUVHI (LOC_Os05g41172), which belongs to the SET 
family and contains an YDG_SRA domain, was found to 
be down regulated. The SRA (SET and RING finger 
associated) domain is believed to play a part in directing 
SUVH proteins to specific chromatin subdomains 
[39-42]. The YDG/SRA domain of KYP/SUVH4 has the 
ability to bind directly to methylated DNA, indicating 
that DNA methylation is necessary for SUVH target- 
ing [43,44]. In Arabidopsis, loss of SUVHI and SUVH4 
causes weak reduction of heterochromatic histone H3K9 
dimethylation [45]. 

In addition, a putative PHD finger protein (LOC_Os07g 
41740) and two RecF/RecN/SMC N terminal domain 
containing proteins (LOC_Os02g04040 (SMC3) and LOC_ 
Osl2g44390 (SMC1)) were up regulated. 

Differential expression of other important proteins 

Many proteins with highly important biological roles were 
also shown to be differentially regulated. The differentially 
expressed proteins included: cleavage and polyadenylation 
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specificity factor, CCAAT/enhancer-binding protein, 
RNA recognition motif containing proteins, OsTOP6B- 
Topoisomerase 6 subunit B protein, DEAD-box ATP- 
dependent RNA helicase, Nucleolar protein NOP5-1, 26S 
proteasome proteins, protease homologue, 14-3-3 proteins, 
importin subunit alpha, DNA topoisomerase 1, cell division 
control protein 48 homolog E, putative Argonaute protein. 

Discussion 

Nuclear proteome and comparison of nuclear protein 
extraction methods 

Proteomic studies on biochemically isolated organelles 
require stringent protein categorization parameters that 
allow for distinction between valid and contaminating 
co-purifying components. In addition, many proteins 
shuttle between the nucleus and cytoplasm and are 
annotated in multiple cellular compartments. There are 
a variety of effective bioinformatic tools for predicting 
nuclear localization based on signal peptides and nuclear 
localization signals, however using these tools for sub 
nuclear domain categorization is not possible. Also, many 
of the entries in the datasets available through these tools 
rely heavily on Uniprot "subcellular localization" field 
keywords. In these circumstances, data available from the 
gene ontology project [28] can be utilized in conjunction 
allowing identified proteins to be classified on their cellu- 
lar localization, biological process, and molecular function. 
Gene ontology is mainly based on available publications, 
which provides relevant evidence of cellular localization. 

Recently, Aki and Yanagisawa (2009) [9] using nanoLC/ 
ESI/MS/MS did extensive studies on the rice nuclear prote- 
ome. Aki and Yanagisawa (2009) [9] possibly identified the 
largest number of nuclear proteins in Rice thus far, using 
co-enrichment with nuclear purification as criteria for 
nuclear localization. In this study, cellular localization was 
classified based on prior publications and GO annotations. 
In recent years, electronic annotation has significantly im- 
proved in terms of specificity, reliability, and coverage [46]. 
Using this cellular localization criterion and two or more 
peptide match for protein identification, we successfully 
identified 382 nuclear proteins. Many of the proteins have 
not been reported in prior nuclear proteome studies. 

We compared the nuclear proteomes extracted by 
phenol and sulfuric acid. The phenol extraction method 
identified 251 nuclear proteins in the nuclei derived 
from protoplasts and 115 proteins in the nuclei derived 
from suspension cells. In contrast, the acid extraction 
identified 137 nuclear proteins in protoplast nuclear 
sample and 165 nuclear proteins in suspension cell nu- 
clear sample. The acid extracted proteins were mainly 
histones, nucleolar proteins, and ribosomal proteins. On 
the other hand, the proteins identified by phenol extrac- 
tion were more diversified. Interestingly, we found that 
further fractionating the phenol extracted proteins by 



sulfuric acid uncovered nuclear proteins that were not 
identified by either method. Sulfuric acid re-extraction 
identified 113 nuclear proteins in protoplast nuclei and 
144 proteins in suspension cell nuclei. Among them, 32 
and 94 proteins were not identified by phenol extraction 
alone of the protoplast and suspension cell nuclei, re- 
spectively. Similarly, 38 and 58 of the proteins were not 
identified in acid extracted protoplast and suspension 
cell samples, respectively. The results suggested that the 
nuclear proteome is highly complex, further fraction- 
ation of the subproteome by acid can lead to a better 
coverage of the nuclear subproteome. Combining phe- 
nol, acid, and their double extraction, we identified 382 
nuclear proteins with two or more peptides, including 
26 transcription factors. The plant (rice) nuclear prote- 
ome has been studied extensively by many authors in 
tissues including rice seedlings, rice suspension cells, 
and rice seed endosperm and evolutionarily conserved 
and glucose responsive nuclear proteins have been iden- 
tified among many other nuclear proteins [9,12,13]. Al- 
though the nuclear purification steps presented all 
appeared to be convincing, the coverage of nuclear pro- 
teins, particularly the low abundant nuclear proteins 
such as transcription factors, remains to be improved. 
Our results suggested that due to the complexity of the 
nuclear subproteome and the presence of high abundant 
proteins such as ribosomal proteins, further fraction- 
ation of the nuclear proteome is necessary to achieve a 
deeper coverage of the nuclear subproteome. 

Regulation of chromatin structure and histone 
modification change in response to cell wall removal 

Previous studies find that removal of the cell wall is con- 
comitant with substantial chromatin reorganization. 
Western blots and isotope-labeling assisted quantitative 
mass spectrometry analyses reveal that the chromatin 
reorganization is associated with substantial histone modi- 
fication changes. Particularly, the H3K18 and H3K23 
acetylation are substantially induced upon removal of the 
cell wall [3]. We are interested in identifying proteins in- 
volved in chromatin reorganization and histone modifica- 
tions. In this study, we found that a histone deacetylase 
(LOC_Os07g06980) was up regulated and a histone-lysine 
N-methyltransferase (LOC_Osllg38900) was down regu- 
lated. Examining the relationship between the regulation 
of these two proteins in response to cell wall removal and 
the histone modification changes caused by cell wall 
removal is of interest [3]. To investigate whether there is 
any causal relationship between the two observed effects, 
we can use the mutants of these differentially regulated 
genes to examine their cellular response to cell wall 
removal and test acetyltransferase activity in response to 
cell wall removal. The structural maintenance of chro- 
mosomal (SMC) proteins function together with other 



Mujahid et al. Proteome Science 2013, 11:26 
http://www.proteomesci.eom/content/1 1 /I /26 



Page 10 of 13 



proteins in a range of chromosomal transactions, includ- 
ing chromosome condensation, sister-chromatid cohesion, 
recombination, DNA repair and epigenetic silencing of 
gene expression [47]. The RecF/RecN proteins are re- 
quired for DNA repair and homologous recombination. 
We found that two RecF/RecN/SMC N terminal domain 
containing proteins structural maintenance of chromo- 
somes (SMC) SMC3 and SMC1, respectively (LOC_OsO 
2g04040 and LOC_Osl2g44390) were up regulated upon 
removal of cell wall. Their potential role in chromatin 
reorganization upon removal of cell wall is worthy of 
further examination. We found that several core histone 
domain containing proteins were up regulated in response 
to cell wall removal. Although the function of this group of 
genes is still unknown, it is a group of very interesting genes 
which should be further explored. A remarkable question is 
whether these proteins are directly involved in the chroma- 
tin re-organization induced by cell wall removal. 

Differentially expressed regulatory proteins and 
cellular process 

To understand the cellular response to cell wall removal 
and the underlying regulatory mechanism, it is essential 
to elucidate the gene regulatory network. Transcription 
factors are the key regulators in gene expression control. 
We found that several transcription factors and tran- 
scriptional regulatory genes are differentially expressed 
in response to cell wall removal. These include two up 
regulated zinc finger proteins and two down regulated 
zinc finger proteins. Other differentially expressed tran- 
scription factors include Helix-loop-helix DNA binding 
protein, factor TF2 (LOC_Os05g03740) containing a 
myb-like family domain, and putative transcription fac- 
tor (LOC_Os09g27850). Our study clearly demonstrated 
differential expression of transcription factors at the pro- 
tein level in response to cell wall removal. In addition, 
we also observed protein level changes in putative DNA- 
directed RNA polymerase and other transcriptional reg- 
ulators or co-regulators. Our results are consistent with 
the dramatic transcriptome change observed in response 
to cell wall removal revealed by oligo microarray studies 
in rice [4]. 

In addition to differential expression of proteins in- 
volved in the transcription process, we also observed 
protein differential expression in RNA binding proteins, 
RNA splicing proteins, ribosomal proteins, translational 
elongation factors, molecular chaperones, protein modi- 
fication proteins, protein degradation proteins. The re- 
sults suggested that the cells responded to cell wall at all 
levels. To further define the regulatory network, we car- 
ried out gene ontology analysis. GO analysis indicates 
that the biological processes tightly associated with cell 
wall removal includes chromatin assembly, nucleosome 
assembly, macromolecular complex subunit organization, 



protein-DNA complex assembly, and DNA packaging. 
Our results clearly indicate that removal of cell wall im- 
poses a tremendous challenge to the cells. Consequently, 
plant cells respond to removal of cell wall in all major 
cellular components and biological processes. 

Materials 

Cell culture 

The rice (Oryza sativa) suspension culture line OC was 
used for all experiments in this study [18,19]. Line OC 
was grown in the dark at 24°C in a gyratory shaker 
under a constant speed of 150 rpm in liquid B5 organic 
medium (pH 5.7) supplemented with 20 g/L sucrose, 
0.5 g/L MES, 2.0 mg/L 2-4-dichlorophenoxyacetic acid 
(2,4-D) as previously reported [3,20,48]. Weekly subculture 
was performed at a dilution of 1:5 (cells:fresh medium). 

Methods 

Protoplast isolation and cell wall regeneration 

OC cells were harvested five days after subculture for 
protoplast isolation. Protoplast isolation was performed 
as previously described [3,20] . Briefly, suspension cells were 
suspended in filter-sterilized enzyme solution containing 
2.5% Cellulase RS (Onozuka RS), 1% Macroenzyme R10 
(Research Products International), 0.4 M mannitol, 80 mM 
CaCl 2 , 0.125 mM MgCl 2 , 0.5 mM MES, and B5 organic 
medium with 2.0 mg/L 2,4-D (pH 5.6). After an incubation 
period in the dark for nine hours at 25°C, the protoplasts 
were collected by first filtering the enzyme solution 
through a 25 (im stainless steel sieve and then centrifuging 
the filtered solution at 120 x g for 5 min. The suspension 
cells were washed several times with protoplast suspension 
medium (0.4 M mannitol, 80 mM CaCl 2 , 0.125 mM 
MgCl 2 , 0.5 mM MES, and 2 g/L N-Z- Amine A in B5 or- 
ganic medium plus 2.0 mg/L 2,4-D at pH 5.6). After proto- 
plasts were washed, they were cultured in sealed petri 
dishes using protoplast suspension medium at a density of 
5 x 10 5 cells/ml in complete darkness at 25°C without 
agitation before being harvested for further study. 

Analysis of new cell wall formation 

New cell wall formation was evaluated by monitoring the 
fluorescence of Fluorescent Brightener 28 (Calcofluor 
White M2R, Fluostain I, Sigma Aldrich, St. Louis, MO) 
using a confocal laser scanning microscope (CLSM) 
Zeiss Axiovert 200 M (Zeiss, Germany) as previously 
described [3]. 

Nuclei isolation and purification 

Rice suspension cells were suspended in nuclear isolation 
buffer (NIB: 10 mM Tris pH 8.0, 2 mM MgCl 2 , 1 mM 
CaCl 2 , 1 mM EDTA, 0.25 M sucrose, 0.1 mM spermidine, 
0.5% Ficoll, 0.5% Triton-X 100 [added freshly], and 1 mM 
PMSF [added freshly], 1 mM DTT [added freshly]). The 
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suspended cells were added to a pre-chilled blender and 
blended on high for 30 seconds. The homogenized slurry 
was first filtered through two layers of cheesecloth, and 
then filtered through a 25 urn stainless steel sieve to 
remove any unbroken cells. The filtered solution was 
centrifuged at 500 x g for 10 min at 4°C. The resulting 
pellet was re-suspended in NIB, under constant shaking at 
4°C for 15 min, followed by centrifugation. Wash steps 
with NIB were repeated three times, followed by layering 
solution on a 2 M sucrose gradient, and centrifugation at 
6000 x g for 30 min at 4°C to pellet purified nuclei. The 
resulting pellet was washed with NIB and used for further 
study. Protoplast nuclei were isolated the same way as 
previously described [3,20]. 

Microscopic observation of purified nuclei 

After purification, the integrity of isolated nuclei was 
assessed by staining with 4', 6' - diamidino-2-phenylindole 
hydrochloride (DAPI). A small volume of the purified 
nuclei was stained with DAPI (0.5 (ig/ml) for 5 minutes 
and images were taken under a D API-filter. 

Nuclear protein extraction 

The protein extraction method is a modification of our 
previous nuclear protein extraction procedure. The pro- 
teins for suspension cell nuclei and protoplast nuclei were 
extracted using phenol extraction as previously described 
[3,49,50]. Three biological replicates were extracted for 
both suspension cell nuclei and protoplast nuclei samples. 
The resulting pellets were further extracted using the acid 
extraction method or directly re-suspended in 8 M urea 
lysis buffer for trypsin digestion. Acid extraction for desig- 
nated nuclear pellets was carried out as previously de- 
scribed [3,51]. To further fractionate the phenol extracted 
proteins, the phenol extracted pellet was suspended in 
0.4 N sulfuric acid and incubated for 2 hours at 4°C with 
constant rotation. After incubation, the solution was 
centrifuged at 16,000 x g for 15 min at 4°C; the resulting 
supernatant was collected and precipitated with a final 
concentration of 33% trichloroacetic acid (TCA) for 
30 min. The TCA precipitated pellet was washed with 
acetone and vacuum dried, followed by suspension in 8 M 
urea lysis buffer. Protein quantification was carried out for 
all samples using the RC DC™ Protein Assay Kit. Three 
replicates were performed for each nuclear protein extrac- 
tion procedure (2 treatments x 3 methods x 3 replicates), 
resulting in a total of 18 mass spectrometric runs. 

Western blot analysis of purified nuclear proteins 

Proteins were separated on a 12% SDS-PAGE gel and 
electrotransfer of gel proteins onto a PVDF membrane 
(Millipore) was carried out at 0.8 mA/cm 2 gel area for 
1.5 hours. Nonspecific binding sites on the membrane 
were blocked overnight with block solution (5% m/V 



non-fat milk, 0.05% v/v tween-20, and 1 X TBS). After 
blocking, the membrane was incubated with respective 
primary antibody for 2 hours at room temperature, 
followed by incubation with respective alkaline phos- 
phatase conjugated secondary antibody for 90 minutes. 
Signal detection was carried out using NBT/BCIP detec- 
tion system. 

RNA Isolation and RT-PCR analysis 

Total RNA was extracted from suspension cell and 
protoplast nuclei using Trizol following manufacturer s 
instructions provided by Invitrogen (Invitrogen, http:// 
www.invitrogen.com). Reverse transcription of RNA was 
performed as previously described with minor modifica- 
tion [52,53]. Rice ubiquitin gene was used as an internal 
control. Resulting PCR products were examined using 
1% agarose gel electrophoresis. PCR primers used in the 
study are supplied in Additional file 4. 

Quantitative real-time PCR Analysis 

Real-time quantitative PCR analysis was performed as pre- 
viously described [53]. The rice ubiquitin gene was used 
as an internal control. The 2" AACT method was used to 
calculate relative transcript levels [54]. Primers used in the 
study are provided in Additional file 4. 

Protein digestion and shotgun proteomic analysis 

Protein digestion was carried out as previously described 
[3,20,25]. Briefly, after dissolving proteins in 8 M Urea 
lysis buffer (pH 7.8), proteins were reduced with 10 mM 
DTT for 1 hour and alkylated with 50 mM IAA for 
1 hour. Subsequently, the urea concentration was re- 
duced to less than 0.6 M for trypsin digestion. Trypsin 
(Promega) was added at a final ratio of 1:50 (protease: 
protein) and digestion was carried out at 37°C overnight. 
Trypsin was inactivated by decreasing the pH to less 
than 2 by adding 2 ul of formic acid. Peptide mixtures 
were desalted with a Michrom Bioresources peptide 
desalting macrotrap following manufacturer s instructions. 
The eluted peptides were vacuum-dried and resuspended 
in 20 ul 5% Acetonitrile, 0.1% formic acid for ID liquid 
chromatography-electrospray ionization tandem MS (ID 
LC ESI MS/MS) using a Surveyor HPLC (Thermo) in-line 
with an ESI ion trap mass spectrometer (LCQ Deca XP 
Plus, ThermoElectron). A reverse-phase column (BioBasic 
C18 column (Thermo 72105-100266)) was used for pep- 
tide separation at a flow rate of 500 nl min-1. Peptides 
were loaded with 5% ACN, 0.1% formic acid for 20 min. 
The elution gradient (all solvents containing 0.1% formic 
acid) was as follows: 5-25% ACN in 450 min, followed by 
25-50% in 130 min, followed by a 20 min wash with 95% 
ACN and then equilibration with 5% ACN for 55 min. 
The extended gradient time was used to compensate for 
the slow scan rate of the instrument. Data was collected 
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over a total duration of 655 min using repetitive MS scans 
directly followed by three tandem MS/MS scans on the 
three most intense precursor masses from the full scan. 
Dynamic mass exclusion windows were 2 minutes long. 
The mass spectra and tandem mass spectra were searched 
against the Oryza sativa non-redundant protein database 
(TIGR, V7.0) downloaded on 1/19/2012 from TIGR Rice 
Genome Annotation (http://rice.plantbiology.msu.edu) by 
using TurboSEQUEST, Bioworks Browser 3.2 (Thermo 
Electron Corp). The database contained 66 338 protein 
entries. Criteria, parameters, and procedure used for pro- 
tein identification were identical to what was previously 
reported [3]. The allowance for missed cleavages was one. 
The peptide (precursor) ion mass tolerance was 1.0 Da, 
and the fragment ion (MS2) tolerance was 0.5 Da. The 
requirement for protein identification was two peptides 
from a protein to meet the following criteria: X- 
correlation >1.9 (+1 charge), >2.2 (+2 charge), >3.75 (+3 
charge); delta correlation value >0.08; probability <0.01. 
Using the reverse database functionality in Bioworks 3.2, 
the peptide and protein false discovery rates were esti- 
mated using the same search criteria as described above 
against the reverse O. sativa database. 

Protein quantification 

TurboSEQUEST (Bioworks Browser 3.2, Thermo Elec- 
tron Corporation), commercial software commonly used 
in mass data analysis was used to generate X corr values. 
The ZX corr quantification method used was as reported 
by Nanduri and Bridges [26,27]. The ProtQuant software 
[27] was downloaded from AgBase [23] database tool 
box (http://www.agbase.msstate.edu/). Quantitative ana- 
lysis criteria and procedure were identical to previously 
reported [3,25]. A peptide X corr value was only considered 
if it passed the following protein identification criteria: 
X-correlation >1.9 (+1 charge), >2.2 (+2 charge), >3.75 (+3 
charge); delta correlation value >0.08; probability <0.01). 
Using the library R statistical package http://www. 
r-project.org/, ProtQuant performed one-way ANOVA 
analysis for proteins identified with three or more peptide 
scans in comparative treatments to determine the statis- 
tical significance of differential expression (p-value). Dif- 
ferential regulation was only considered for proteins with 
a p-value < 0.05. 

Gene ontology annotation 

In order to carry out protein functional categorization, the 
gene ontology (GO) rules provided with the GO browser 
at http://www.geneontology.org/ [28] were followed. Gene 
ontologies can be classified into three independent groups: 
biological process (BP), molecular function (MF), and 
cellular component (CC). Using the GORetriever tool 
available at AgBase [23] (http://www.agbase.msstate.edu/), 
GO annotations were assigned. If GO annotations could 



not be retrieved using this tool, other websites including 
Uniprot, TIGR, NCBI, and Gramene were used to retrieve 
annotations. GOSlimViewer (available at AgBase) tool was 
used to retrieve GoSlim ids. Functional categorization of 
genes was also carried out according to the GO rules [28] 
atagriGO [29,30]. 
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